function [az, ele] = azel_compute(svxyzmat,usrxyz)
% edn下仰角和方位角
ele = [];
az = [];
wgs84 = wgs84Ellipsoid('meter');
% for j = 1:length(svxyzmat)    
    r_lla = ecef2lla(usrxyz);
    [s_n,s_e,s_d] = ecef2ned(svxyzmat(:,1),svxyzmat(:,2),svxyzmat(:,3),...
    r_lla(1),r_lla(2),r_lla(3),wgs84);
    R = sqrt(s_n.^2+s_e.^2+s_d.^2);
    ele = abs(asind(s_d./R));
    az = atan2d(s_e,s_n);
% end